% Test the function makeppft filter.
%
% Yoel Shkolnisky, December 2010.
precond=0;
for n=2:2:32
    
    tic;
        H1=makePPFTFilter3D(n);
    t1=toc;
    
    tic;
    filename=sprintf('filters\\H3Dfilter%dprecond%d.dat',n,precond);
    if exist(filename,'file')
        H2=load3DHFilter(filename);
    else
        H2=save3DHFilter(filename,n,precond);
    end
    t2 = toc;

    H3=H1.filter(2:2*n,2:2*n,2:2*n);
    
    reportTestResult(n,(norm(H2.filter(:)-H3(:))./norm(H2.filter(:))),t1/t2);    
end

